Tailoring of phononic band structures in colloidal crystals 
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We report an experimental study of the elastic properties of a two-dimensional (2D) colloidal crys- 
tal subjected to light- induced substrate potentials. In agreement with recent theoretical predictions 
1 the phonon band structure of such systems can be tuned depending on the symmetry and depth 
of the substrate potential. Calculations with binary crystals suggest that phononic band engineering 
can be also performed by variations of the pair potential and thus opens novel perspectives for the 
fabrication of phononic crystals with band gaps tunable by external fields. 
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Materials with periodic variations in their elastic prop- 
erties receive currently much interest as phononic crys- 
tals. Analogue to light propagation in photonic crystals 
[2j [3] , the transmitted spectrum of sound waves travel- 
ing through phononic crystals exhibit band gaps whose 
frequency is determined by the length scale on which the 
elastic properties are modulated. Fabrication of phononic 
crystals is achieved by embedding regular arrays of elas- 
tic inclusions in an appropriate matrix. Experiments 
with millimeter-sized inclusions [H [5] or sub-micron-sized 
holes immersed in elastic host materials [6] [7] indeed 
show acoustic band gaps at ultrasonic, sonic and hyper- 
sonic frequencies. Recent Brillouin spectroscopy mea- 
surements on crystals made of sub-micron-sized colloidal 
particles immersed in a liquid matrix demonstrated band 
gaps in the hypersonic regime with the possibility to tune 
the band gap by exchange of the surrounding liquid [8]. 

In this letter, we experimentally investigate the 
phononic properties of a two-dimensional (2D) crystal of 
colloidal particles being subjected to a periodic substrate 
potential. Depending on the substrate strength and par- 
ticle interactions, the phononic band structure and thus 
the position and width of phononic band gaps can be 
largely tuned. Because this concept applies not only to 
micron-sized colloids but also to much smaller particles, 
this suggests tailoring the phononic properties of atoms 
or molecules confined to extended optical lattices [9]. 

Experiments were performed with an aqueous suspen- 
sion of highly charged polystyrene spheres with diameter 
<j = 2Apm and a polydispersity below 4%. The parti- 
cles interact via a screened Coulomb potential <3>(r) ex 
Z 2 exp(— K,r)/r with Z « 10000 the renormalized surface 
charge and n~ l ~ 300nm the screening length. Both val- 
ues were determined according to a procedure described 
in [10]. As sample cell we used a cuvette made of fused 
silica with 200pm spacing between top and bottom plate 
which was connected to a standard closed deionization 
circuit to maintain stable ionic conditions during the 
measurements [11]. After sedimentation, the particles 
form a 2D colloidal system close to the bottom plate. 

One- and two-dimensional substrate potentials were 
created by superimposing two perpendicularly aligned 



one-dimensional periodic interference patterns created 
with a P = bW frequency-doubled Nd:YVC>4 laser 
(wave length = 532nm). The polarizations of the inter- 
ference patterns were adjusted perpendicular, therefore 
they act as two independent ID periodic substrate po- 
tentials for the colloidal particles Ui(x) = Uq cos(2ttx /di) 
with Uq the potential amplitude, di the lattice constant 
and i = 1,2 [12 . Because the potential amplitude Uq 
scales linearly with the intensities of the laser beams [13] 
this allows to continuously adjust the strength of the un- 
derlying laser potentials. An additional laser beam was 
scanned around the central region of the sample to cre- 
ate a boundary box whose size could be continuously ad- 
justed by a pair of computer-controlled galvanostatically- 
driven mirrors. This allowed us to adjust the particle 
density p with an accuracy of Ap/p w 0.01 [T4] . 

First, we continuously decreased the mean particle dis- 
tance a = yj 2/\[3p to approximately 4pm where under 
our salt conditions spontaneous crystallization of the col- 
loidal monolayer occurs. Next, the lattice constants d\ 
and were chosen to meet commensurate conditions, 
i.e. d\ = (v // 3/2)a « 3.5pm and di = a/2 « 2pm (see 
Fig. [I]). Particle positions were determined for different 
combinations of interference pattern intensities {Uq , Uq } 
from sequences consisting of several thousand images us- 
ing digital video microscopy at an acquisition rate of 2 
frames per second [10] . From these data we finally ob- 
tained the particle trajectories r^{t) = (x^(t),y^(t)) with 
p = 1 ... N and t the time. To avoid boundary effects, 
we only considered the central region of the field of view 
(300pm x 200pm) containing N « 1000 particles. 

For the experimental determination of the phononic 
band structure we followed the approach already intro- 
duced in [TJ [I5j [16] which is based on the analysis of 
the dynamical matrix. Briefly, the branch X s (q) of the 
phononic band structure with q the wave vector and po- 
larization s is given by the eigenvalue of the dynamical 
matrix denoted as D a p(q) (a, (3 = x,y). The latter is ob- 
tained from the measured particle displacements u(R, t) 
relative to their lattice sites R 

D a0 {q) = k B T/(u* a (q,t)up(q,t)) t (1) 
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with u(q,t) the Fourier transform of u(R,i) and (. . .) t 
the temporal average. Using the equipartition theorem, 
this expression has been derived within the harmonic ap- 
proximation of the potential energy 



H=1/2J2 <(q)D a p(q)u (q) 
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which has been experimentally confirmed to be valid 
for colloidal systems [15 j. Note, that due to the ab- 
sence of true long-range order in 2D systems, the lat- 
tice sites R have to be determined through temporal av- 
eraging [TTJ [18]. Within the harmonic approximation 
the pair interaction between particles can be modeled 
as springs with the spring constant ko given by the sec- 
ond derivative of the pair potential at the mean par- 
ticle distance ko = [d 2 ®(r)/dr 2 ] r=a . With the corre- 
sponding values Z and k taken from above we obtain 
k « 150k B T/a 2 « 1 • l(T 7 J/ra 2 . The periodic sub- 
strates Uq introduce an additional set of springs which 
pin the particles to their lattice sites R with spring 
constants k{ = UqCZtt / di) 2 . With the laser intensities 
used in our experiments we achieved maximum values of 
ki ~ 2.5&Q- Since the phononic band structure depends 
only on the ratio between pair and particle-substrate in- 
teraction, in the following it will be expressed in units 
of ko. We confirmed the harmonic approximation to 
be valid over the entire range of observed displacements 
(\ u a\ < lfJ>m) by determining the effective single-particle 
potential U e f f(u a ) /ksT = —\n[P(u a )] with P(u a ) the 
normalized particle-displacement distribution. 

We will start our discussion by first describing the re- 
sults obtained when only a single ID substrate potential 
was applied (Uq = 0) (see Fig. [2J. A typical snapshot of 
the colloidal crystal is shown in Fig. [2|a) (the ID laser 
potential which is aligned vertically is invisible because 
it is blocked with optical filters). Fig. [2^b) illustrates the 
irreducible path along which the phonon band structure 
has been analyzed. It corresponds to that route which 
circumvents that part of the first Brillouin zone (gray 
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Figure 1: Illustration of the colloidal system subjected to com- 
mensurate optical interference patterns. Gray dots represent 
colloidal particles and gray lines indicate the minima of the 
light-induced ID periodic potentials. The lattice constants d\ 
and c?2 of the ID periodic potentials are chosen to be com- 
mensurate with the mean particle distance a. 



area) being the smallest repeat unit of the band struc- 
ture. Due to the two-fold symmetry of the ID periodic 
substrate potential, this unit comprises a quarter of the 
first Brillouin zone. Figures [2jc)-(e) and (f)-(h) show the 
2D pair correlation function g(x,y) and the correspond- 
ing phonon band structure for increasing laser intensities 
fci/fco = (spontaneous crystal) (c,f), k\/ko = 1 (d,g) 
and fci/fco = 2.5 (e,h), respectively. Since the particles' 
motion becomes more confined in the direction perpen- 
dicular to the potential troughs, the fluctuations in x- and 
y-direction become asymmetric with increasing (k\/ko) 
as this is observed in the corresponding g(x : ?/)-plots. 




Figure 2: (color online) Experimental realization of a 2D col- 
loidal crystal on a ID periodic substrate [see Fig.[j]a)]. a) Mi- 
crograph of the 2D colloidal crystal in the central region of the 
field of view (300/ira x 200/xra). b) Zero and first order Bragg 
peaks (black dots) and first Brillouin zone (gray area). The 
black solid line represents the irreducible path along which 
the phonon band structure is plotted and the path sections 
are marked by numbers. c)-e) 2D correlation function g(x,y) 
and f)-h) phonon band structure X s (q) along the irreducible 
path for c),f) the spontaneous crystal, i.e. ki/ko = 0, d),g) 
ki/ko — 1, and e),h) ki/ko = 2.5. Closed and open sym- 
bols refer to the upper and lower band, respectively. Solid 
lines represent band structures calculated in the framework 
of harmonic lattice dynamics pQ . 
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From the above it is clear that changes in g(x,y) 
also affect the dynamical matrix and thus the phononic 
band structure. As we will demonstrate in the follow- 
ing, those changes strongly depend on the polarization 
of the phonons. Usually, the polarization s refers to lon- 
gitudinal and transversal modes in the direction of high 
crystal symmetry. Because we also determined the band 
structure in directions with lower symmetry, in the fol- 
lowing we will refer to the different polarizations by their 
vertical position in Figs. |2^f)-(h) as lower (s = I) and 
upper (s = u) bands instead. Interestingly, both bands 
are affected by the substrate potential in a rather dif- 
ferent way. We observe a pure shift of the upper band 
X u {q) (closed symbols) along the sections 1 and 2 while 
the lower band \\{q) (open symbols) is not affected. This 
can be easily understood because X u (q) and \i(q) have 
polarization vectors in the x- and the ^-direction along 
these sections, respectively. As a consequence, X u (q) is 
increased by the spring constant fci/fco of the periodic 
potential. The bands have the same polarization along 
section 4 of the irreducible path. Therefore, the upper 
band is not influenced by the periodic potential and the 
lower band is shifted by fci/fco • We observe a different 
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Figure 3: (color online) Experimental realization of a 2D col- 
loidal crystal in the presence of two ID periodic substrates 
rotated by 90° (see Fig. [I]). The graphs are 2D contour 
plots of the phonon band structure in the first Brillouin 
zone derived from experimental (left column) and analyti- 
cal data (right column), a) {fci/fco, fo/A:o} = {0,0}. b) 
{fci/fco, fe/M = {1.4,0.2}. c) {fci/fco, fc 2 /fco} = {1.5,1.5}. 
d) {fci/fco, fc2/fco} = {0.5,1.5}. Experimental contour plots 
were determined for approximately 800 values. Calculated 
plots were obtained in the framwork of harmonic lattice dy- 
namics using approximately 3000 values pp. 



behavior along section 3; here, both the upper and the 
lower band are influenced because the respective polar- 
ization does not point in the x- and the ^-direction. The 
upper band becomes largely deformed and is almost en- 
tirely flattened in Fig. |2|h). Because the group velocity 
v(q) oc Vy/X(q) vanishes in those regions, this opens in- 
teresting perspectives regarding the tailoring of acoustic 
and thermal properties of 2D crystals. We also compared 
the experimentally determined band structures with cal- 
culations based on the harmonic approximation [1 using 
fci/fco as a fit parameter, the latter being in excellent 
agreement with the corresponding values independently 
obtained from n~ x and Uq. The calculated band struc- 
tures are plotted as solid lines in Figs. |2jf)-(h). 

Similar experiments have been also performed for 2D 
substrate potentials ({^0,^0} 7^ {0?0}) an d demon- 
strate that both the upper and the lower band can be 
individually shifted and deformed by the presence of the 
substrate potentials. Rather than showing a specific path 
through the Brillouin zone, here we have chosen a 2D 
representation of the phononic band structure which is 
shown in Fig. [3] for different combinations of substrate 
strengths. For comparison, we also show calculated 2D 
phonon band structures obtained within the harmonic 
approximation using fci/fco and fc2/&o as fit parameters 
Pp. For the substrate- free case where {fci/fco, fe/^o} = 
{0, 0}, as expected we observe the six- fold symmetry of a 
spontaneous crystal [Fig. [3ja)]. In Fig. [3^b), the spring 
constants are {fci/fco, &2/&0} = {1-4, 0.2}. Here, the sym- 
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Figure 4: Phonon band structure (solid symbols) of a para- 
magnetic 2D binary crystal (see inset) calculated in the har- 
monic approximation. The black particles in the inset have a 
magnetic susceptibility 20 times larger than the susceptibility 
of the gray particles. The magnetic field B is perpendicular to 
the system plane. The band gap can be closed for B = and 
A > if all the particles additionally interact via a Lennard- 
Jones potential (open symbols). 
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metry breaks down to a two-fold symmetry due to the 
dominating k\. Fig. J3^c) corresponds to the situation 
where both interference patterns have identical strengths 
{fci/fco? = {1-5, 1.5} and as a consequence the six- 

fold symmetry is restored. Finally, Fig. [3^d) shows the 
case where {fci/fco, &2/&0} = {0-5, 1.5}. Similar as in Fig. 
(3^b) we find a two-fold symmetry but here rotated by 90° 
because of > fci/fco- 

After having discussed the possiblity to tailor phononic 
band structures by subjecting monolayers to substrate 
potentials, finally, we want to discuss whether phononic 
band engineering can also be achieved by adjustable 
anisotropic pair potentials. Here, we exemplarily con- 
sider a system with magnetic dipole-dipole interactions 
^mfj) = cfi /\r\ 3 — 3(d • r) 2 /|r| 5 . For paramagnetic par- 
ticles the dipole moment d scales as d = xB with x the 
magnetic susceptibility and B the magnetic field. The 
pair interaction becomes anisotropic when the magnetic 
field is not applied perpendicular to the plane of the 
2D crystal. Our calculations based on the approach in 
PQ indeed show that tuning of phononic bands is pos- 
sible under such conditions, however, we also find that 
anisotropic pair potentials are not sufficient to create 
band gaps in monodisperse systems. This is because only 
acoustic branches exist in such systems where X(q) — > 
for q — » 0. To overcome this limitation, one can create 
additional optical branches (X(q) > for q = 0). This 
is achieved by extending the unit cell of the crystal by 
adding particles with different magnetic properties to the 
system (black and gray spheres in the inset of Fig. [4|. 
For our calculations we assumed a magnetic susceptibil- 
ity ratio X1/X2 — 20 between the black and gray particles 
and the magnetic field B orientated perpendicular to the 
system. The calculated phonon band structure is shown 
in Fig. [4] as closed symbols and indeed shows four addi- 
tional optical branches. Most importantly, we observe a 
full band gap whose frequency uj oc y/X and width scale 
linearly with the magnetic field. As a consequence, the 
band gap can only be closed in the limit B — > where 
both the pair interaction and the corresponding spring 
constants X(q) vanish; therefore, we considered an addi- 
tional pair interaction which was exemplarily assumed as 
a Lennard- Jones potential </>z,j(r) oc [(cr/r) 12 — 2(<j/r) 6 ]. 
Similar as above the band gap can be completely closed 
at B = but now at finite frequencies, here determined 
by the Lennard- Jones interaction (see open symbols in 
Fig. [4|. This would allow to extend the possibility to 
tailor phononic band structure also to situations where 
surface potentials can not be modified. Preliminary cal- 
culations indicate that a similar behavior can be observed 
in three-dimensional systems. 

In conclusion, we have experimentally demonstrated 
that the phonon band structure of a 2D colloidal crystal 
can be greatly tuned by the strength of periodic substrate 
potentials. Depending on the symmetry of the applied 



substrate potential different phonon polarizations can be 
tailored rather independently. Calculations with binary 
crystals of paramagnetic particles with adjustable pair 
interactions indicate that phononic band engineering can 
be also performed in situations where substrate poten- 
tials can not be tuned and thus opens novel perspectives 
for the fabrication of phononic crystals with band gaps 
adjustable by external fields. 
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